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Evidence for a 2:3 resonance in Sco X-l kHz QPOs 
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Abstract. We find evidence that the two high frequency QPOs in Sco X-l are, more often than not, approximately 
in the 2:3 frequency ratio familiar from studies of black hole candidates (e.g., XTE J1550-564, Remillard et al. 
2002). This implies that the double kHz QPO phenomenon in neutron stars has its origin in properties of strong- 
field gravity and has little to do with the rotation of a stellar surface or any magnetic field structure anchored in 
the star. 
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1. The larger picture 

Recent X-ray observations showed that in some black 
holes, as in neutron stars, high frequency quasi-periodic 
oscillations (QPOs) come in pairs. The frequencies of these 
pairs in the first three black hole sources where they have 
been observed are 300 Hz and 450 Hz in GRO J1655-40, 42 
Hz and 70 Hz in GRS 1915+105, and 184 Hz and 276 Hz 
in XTE J1550-564 (Strohmayer 2001a,b; Remillard et al. 
2002). Clearly, they are in a 2:3 ratio in two of the sources, 
and in a 3:5 ratio in the third (Abramowicz and Kluzniak 
2001; Kluzniak and Abramowicz 2001,2002; Remillard et 
al. 2002). In general relativity, characteristic frequencies 
scale inversely with the mass. If one scales the 42 Hz 
QPO frequency of the ~ 20M Q black hole candidate GRS 
1915+105 to about one solar mass, i.e., the mass of a neu- 
tron star, one obtains ~ 800 Hz, a frequency typical of 
the lower of the two "kHz" QPO frequencies detected in 
neutron stars. This in itself suggests that the QPO phe- 
nomenon in black holes and in neutron stars may have the 
same origin. We think that the frequency ratios may be 
even more revealing. 

In this paper we examine the published frequency data 
of the well-studied source Sco X-l, a prototype of bright 
neutron-star candidates, to see whether the two "kHz" 
QPO frequencies may be in the ratio of two low integers, 
such as 2:3 or 3:5, as they are in black hole systems. If 
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such a similarity were well established between neutron 
stars and black holes, it would imply that the "twin" high 
frequency QPOs originate in a mechanism operating in- 
dependently of the presence of any stellar surface with its 
co-rotating magnetic fields. If the high frequency QPOs 
are fundamentally similar in neutron stars and in black 
holes, they are all likely to be a manifestation of strong- 
field gravity. 



2. The inner disc radius in neutron-star LMXBs 

Is the accretion process in neutron star members of low- 
mass X-ray binaries (LMXBs) radically different from that 
in their black hole counterparts, or are the two processes 
fundamentally similar? This question predates the discov- 
ery of high frequencies in these systems. One school of 
thought held that neutron stars in LMXBs have magnetic 
fields and spin rates comparable to that of millisecond ra- 
dio pulsars, which would imply that the accretion disk is 
terminated far from the stellar surface, and that the high- 
est orbital frequencies would be at most two or three hun- 
dred Hertz (within 50 Hz or so of the presumed stellar spin 
rate). Others suggested that the magnetic field in these 
systems may be so low as to play no dynamical role, and 
that the discs (presumed to be geometrically thin) may be 
terminated close to the central compact object by effects 
of strong gravity much as they are in black holes. On this 
(minority) view, orbital frequencies in the kilohertz range 
were expected (Kluzniak, Michelson, and Wagoner 1990). 
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It would appear that the discovery of kHz quasiperiodic 
oscillations in neutron star systems favors this latter view. 

The kHz QPOs have fairly similar properties in a wide 
range of systems, in particular the highest frequencies have 
similar values in X-ray bursters and the much brighter 
sources, such as Sco X-l (see van der Klis, 2000, for a 
review). This argues against any special mechanism re- 
lated to accidental properties of the system, such as the 
value of the magnetic field or of the accretion rate (al- 
though many models based on such finely tuned proper- 
ties have been proposed). If, then, the oscillation arises in 
the accretion flow in the disk, the orbital frequency sets 
the scale of expected oscillation frequencies (Bath 1973). 
In fact, the relevant frequency is the meridional epicyclic 
frequency It is identical to the Keplerian frequency for 
Newtonian 1/r gravity, but for accretion disks around 
neutron stars in general relativity the meridional epicyclic 
frequency is typically somewhat smaller than the orbital 
frequency. Observations of a kHz frequency then strongly 
disfavor a geometry of accretion in which the maximum 
orbital frequency is only a few hundred Hertz. This seems 
to rule out surface magnetic dipole fields in excess of about 
10 8 G (Kluzniak 1998), and tends to favour the view that 
the accretion disk in these neutron star systems termi- 
nates close to the marginally stable orbit of general rela- 
tivity (Lipunov and Postnov 1984; Kluzniak and Wagoner 
1985; Sunyaev and Shakura 1986; Kluzniak, Michelson, 
and Wagoner 1990; Kluzniak and Wilson 1991; Biehle 
and Blandford 1993; Zhang, Strohmayer, and Swank 1997; 
Kaaret, Ford and Chen 1997; Kluzniak 1998; Zhang et al. 
1998). 

3. Is the stellar spin frequency directly observed 
in kHz QPOs? 

Patterson (1979) pointed out that if there is direct inter- 
action between a structure corotating with the star and 
a Keplerian accretion disk, a QPO at the beat frequency 
|1/P* — 1/-Fk| may result. Here P* is the period of stellar 
rotation and Pk the period of orbital motion. Patterson 
was discussing the white dwarf system AE Aqu, and the 
co-rotating structure was an "illuminating searchlight." 

After the discovery of ~ 40 Hz QPOs in LMXBs, this 
idea was adapted to a neutron star rotating at a presumed 
frequency of ~ 160 Hz and beating close to a presumed 
Alfvcn radius against a ~ 200 Hz Keplerian frequency at 
the inner edge of an accretion disk (Alpar and Shaham 
1985), and given the name "beat-frequency model" (Lamb 
et al. 1985). But a maximum 200 Hz Keplerian frequency 
is incompatible with the more recently discovered kHz os- 
cillations in the same sources, so the name has been duly 
changed to "sonic-point model" with one QPO still cor- 
responding to a (greatly increased) beat frequency and 
the other peak to the orbital frequency (Miller, Lamb and 
Psaltis 1998). It has been noted (e.g., Mendez et al. 1998; 
Jonker et al. 2002a) that, in fact, the difference between 
the two observed high frequencies in neutron star LMXBs 
is not constant, as it presumably should be in the beat fre- 



quency model. Accordingly, other models were proposed 
in which the interaction between a Keplerian disk and the 
stellar magnetosphere ("co-rotating structure") gave rise 
to frequencies in greater concordance with the observa- 
tions (Osherovich and Titarchuk 1999). 

Alas! The discovery of double high frequency QPOs 
also in black hole candidates (Strohmayer 2001a, b; Miller 
et al. 2001; Remillard et al. 2002) has cast a deep shadow 
over the notion that there may be a direct relationship be- 
tween stellar spin and the presence of two high frequency 
QPOs. After all, black holes do not have a surface and 
their spin frequency cannot be directly observed. 

The general belief seems to be that the distribution of 
frequency ratios is fairly uniform, as it is for the source 
GX5-1 (Fig. 1). Here, we note that this is indeed the ex- 
pectation in any model, where one frequency varies freely 
over some range, and the second differs from the first by 
a nearly constant value. There would be no reason for the 
observed frequency ratio to cluster about any particular 
value, much less a value equal to the ratio of two integers — 
in the "co-rotating structure" model the ratio between the 
constant stellar rotation rate and the orbital frequency is 
accidental and cannot be the ratio of two small integers. 
Even stranger would be a clustering about a value close to 
that observed in black hole systems, e.g., a 2:3 ratio. And 
yet, such a clustering is exactly what we find in Sco X-l 
(Fig. 2). 

One weakness of our argument is the uncertain nature 
of Sco X-l. Not being an X-ray burster, Sco X-l enjoys 
the status of a neutron star candidate. We have performed 
a similar analysis of other neutron star sources, and the 
peak at about 2/3 ratio is present in those data as well 
(Fig. 3). 

4. Evidence for the 2:3 frequency ratio in Sco X-l 

The high frequency QPOs of Sco X-l have been discussed 
in detail by van der Klis et al. 1997. In particular, the 
frequencies of the two QPOs have been presented for a 
number of separate continuous observations of 1 to 3 ks 
each. The frequency, f 2l of the upper QPO ranged from 
about 850 Hz to slightly less than 1100 Hz. The other 
QPO had a frequency, /i, lower by about 300 Hz to 200 
Hz, respectively. As usual, the difference was clearly 

anticorrelated with the frequency /2. 

We have computed the ratio /1//2 for each observation 
of Sco X-l reported in van der Klis et al. 1997. We find it 
to be in a much narrower range than that of the combined 
data of all neutron star sources (Fig. 3). 

For the Sco X-l data considered, 0.64 < /i// 2 < 0.8. 
We have prepared a histogram of the distribution of this 
QPO frequency ratio, by binning the ratio in suitable in- 
tervals (Fig. 2). It is at once evident that the distribution 
is not uniform, it is instead sharply peaked about a central 
value fairly close to a 2:3 ratio (indicated in the figure by 
a vertical dashed line). This conclusion is not affected by 
the bin width, and is fairly representative of the complete 
sample. 
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We are not aware of any observational bias which 
would favor a selection of those data in which the QPO 
frequency ratio would cluster about 0.67, neither are the 
observers (P. Jonkcr, private communication). 

5. Statistical significance of the peaks in the 
frequency ratio 

To quantify our claim of clustering about a preferred value 
of the frequency ratio of "kHz" QPOs in Sco X-l, we have 
fit to the observed distribution of frequency ratios a suit- 
ably normalized Lorentzian as a trial probability distribu- 
tion, 

dp _ N(r ,X) 
dr (r — r G ) 2 + A 2 ' 

where ro and A are free parameters, and N(r ,X) is the 
normalization factor. Here, r = fi/f%. We use the max- 
imum likelihood method (Bulik, Kluzniak, and Zhang 
1999) to find the probability distribution that best de- 
scribes the data, i.e., the most likely value of the two 
parameters. Once these are determined, we calculate the 
probability that the data has been drawn from this distri- 
bution using the Kolmogorov Smirnov test. The likelihood 
function is the product of probabilities of measuring each 
individual frequency ratio £(ro,A) = Hi(dp/dr\ ri ). 

We searched the two dimensional parameter space for 
the best values of ro and A and found a strong peak 
at ro = 0.68 and A = 0.018. The Kolmogorov Smirnov 
test shows that the probability that the data has been 
drawn from the distribution described by eq. (1) with 
these two values of the parameters is p = 0.39. This con- 
firms that our initial choice of the shape of the distribution 
(Lorentzian) is close to the original distribution of the ra- 
tios, and that the probability distribution maximizing the 
likelihood function is acceptable. 

In Figure 2, we superpose on the histogram this best- 
fit Lorentzian distribution (smooth curve), normalized to 
the total number of events. The offset of the peak from 
the value 2/3 is statistically significant. A fit with a forced 
value of the centroid at ro = 2/3 yields a slightly broader 
distribution with A = 0.021, but the K-S probability is 
unacceptable, at p = 7.1 x 10 . 

The data of fig. 2, here binned in intervals of width 
dr = 0.02, exhibit an additional peak around the value 
of r ks 0.78. We have also tried fitting two Lorentzians. 
The best fit values are r = 0.679 and 0.7763 for the cen- 
troids, and A = 0.01369 and 0.0049 for the widths of the 
two peaks, respectively. Although inclusion of a second 
Lorentzian improves the fit (K-S probability of 0.801), 
from the statistical point of view the second peak is not 
strongly required by the data, as the K-S test probability 
is already acceptable with just one Lorentzian. 

A careful analysis of Fig. 3 in van der Klis et al. 1997, 
especially panel d), suggests that the second peak may be 
real — the power ratio of the two QPOs departs from the 
general trend when /1//2 becomes greater than 0.75 (i.e., 
when f 2 > 1050Hz). From the data (Fig. 2 of this work 
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Fig. 1. The observed distribution of frequency ratio of the 
two kHz QPOs of GX5-1. Also shown is the location of 
the 2:3 ratio (dashed line). The data is from Jonker et al. 
2002b. 




Fig. 2. The same for Sco X-l. The data is from van der 
Klis et al. 1997. Also shown is the best fitting Lorentzian 
(eq. [1]), note the slight offset of the peak from the value 
2/3. 

and Fig. 3 of van der Klis et al. 1997), it would appear 
that when the QPO frequency ratio is /1//2 = 2/3 (to 
within a few percent), the upper QPO has about twice as 
much power as the lower QPO, and the upper QPO power 
steadily decreases with increasing frequency ratio /1//2, 
until /1//2 ~ 0.75, when the upper QPO again becomes 
twice, and even three times, as powerful as the lower QPO. 

Recent observations of the burster 4U 1636-53, support 
the reality of both peaks in Fig. 2. The frequency differ- 
ence between the QPOs in 4U 1636-53 seems to undergo an 
abrupt transition from a peak separation of about 330 Hz 
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Fig. 3. The combined distribution of frequency ratios 
for 4U 0614+09, 4U 1608-52, 4U 1636-53, 4U 1702-43, 
4U 1705-34, 4U 1728-34, 4U 1735-44, 4U 1820-30, 4U 1915- 
05, Cyg X-2, GX 17+2, GX 340+0, GX 349+2 (Sco X-2), 
KS 1731-26, and XTE J2123-058. 



VL Z , in that orbit in which the radial epicyclic frequency 
is fi r = (2/3)fi z . For the Schwarzschild metric this oc- 
curs at r2:3 = 5.4r g = 16.2km/MQ, and comparison with 
the observed frequency of 900 Hz in Sco X-l implies that 
M — IMq, with modestly higher values for rotating neu- 
tron stars. 

The origin of the second peak of frequency ratios, 
at ~ 7/9, remains a puzzle, as does the flat distribu- 
tion of frequency ratios seen in GX5-1 and some other 
sources. Detailed comparison of observations and the the- 
ory of parametric epicyclic resonance for realistic space- 
time metrics of rapidly rotating neutron stars will be the 
subject of another paper. 
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On on^ side of the transition the frequency ratios are 
0.662, 0.680, 0.685, 0.725, i.e., they may well have been 
sampled from the 2/3 peak in the distribution of Sco X-l; 
on the other side, the ratios vary between 0.75 and 0.78. 

6. Interpretation of the results 

Figs. 1, 2, and 3 make it clear that, unlike GX5-1, most 
neutron star sources do not show a flat distribution of kHz 
QPO frequency ratios. It is hard to escape the conclusion 
that a resonance is responsible for the distribution of fre- 
quency ratios in Sco X-l. The 2/3 value, identical to the 
ratio seen in the black hole candidates GRO J1655-40 and 
XTE J 1550-564 (whose measured masses are ~ 6M Q and 
~ SMq) is a very strong indication that the high frequency 
QPOs are a manifestation of strong-field gravity, and not 
of the presence of a stellar surface with any co-rotating 
structure. 

Our argument about the role played by strong field- 
gravity in forming the high frequency QPOs is based on 
the similarity between neutron stars and black holes, for 
this we need to assume that the microquasars GRO J1655- 
40 and XTE J1550-564 are in fact accreting black holes. 

To make the argument secure, one needs a theory of 
the origin of high frequency QPOs able to accomodate 
both the sharp frequency ratios in the observed black 
holes and the broader distribution observed in Sco X-l 
and some other neutron star sources. It must explain a 
preference for a 2/3 ratio of the high frequency QPOs 
and the offset between the peak of the distribution of 
Fig. 2 and the value 2/3. We note that all these facts seem 
to be in accord with the theory of epicyclic parametric 
resonance in relativistic accretion disks. On that theory 
(Kluzniak and Abramowicz 2002), the higher frequency 
QPO occurs at the frequency of meridional oscillations, 
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